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Abstract 

Using a coherent state representation we derive many-body probability distributions and 
wavefunctions for the Chern-Simons matrix model proposed by Polychronakos and compare 
them to the Laughlin ones. We analyze two different coherent state representations, cor- 
responding to different choices for electron coordinate bases. In both cases we find that 
the resulting probability distributions do not quite agree with the Laughlin ones. There is 
agreement on the long distance behavior, but the short distance behavior is different. 
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1. Introduction 



There has recently appeared an interesting connection between quantum Hall effect and 
noncommutative field theory. In particular Susskind proposed in , that the Laughlin states 
at filling fractions v = for a system of an infinite number of electrons confined in the 
lowest Landau level (LLL) can be described by a noncommutative U{1) Chern-Simons the- 
ory. The fields of this theory are infinite matrices which act on an infinite Hilbert space, 
appropriate to account for an infinite number of electrons. In the same spirit, Polychron- 
akos, later, proposed a finite matrix model 0, a regularized version of this noncommutative 
Chern-Simons theory in an effort to describe finite systems, of limited spatial extent with 
a finite number of electrons. This matrix model was shown to reproduce the basic features 
of the quantum Hall droplets and their corresponding excitations , such as boundary and 
quasihole excitations at filling fraction u = 2f^- 

In a subsequent paper 0, Hellerman and Raamsdonk, trying to make the connection 
between quantum Hall effect and the noncommutative matrix model more transparent, ana- 
lyzed the states of the theory and concluded that the states of the matrix model are in one- 
to-one correspondence with the Laughlin states describing QHE at filling fraction u = 2ffi- 
Similar arguments were put forward for the excited states. Although this is an interesting 
observation, the existence of a one-to-one mapping between states is not enough to prove 
the equivalence of the two theories. 

Since the mapping in is somewhat formal at the level of states, in this paper we try 
to go one step further and compare the two theories at the level of the wavefunctions. This 
requires a notion of coordinates, which is introduced via a coherent state representation for 
the states of the matrix model. We will present two different ways of arriving at a coherent 
state representation, which are not equivalent in general. In the special case of u = 1 { 
^ = 0), it turns out that both coherent state representations of the corresponding matrix 
model reproduce the v = 1 Laughlin wavefunction. For z/ = though, the two ways 
produce different probability distributions and neither agrees with the Laughlin one. The 
probability distributions emerging from the different coherent state representations have a 
long distance behavior similar to the corresponding Laughlin one, but quite a different short 
distance behavior. 

This paper is organized as follows. In section 2 we give a brief discussion of fermions in 
the lowest Landau level and the Laughlin wavefunctions. In section 3 we briefly review the 
Chern-Simons flnite matrix model proposed by Polychronakos. In section 4 we present two 
different approaches towards a coherent state representation of the matrix states. In section 
5 we comment on our results. 

2. Fermions in the lowest Landau level 

It is well known by now that the two-dimensional configuration space of charged particles 
in the LLL is equivalent to a one-dimensional phase space and therefore noncommutative 
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0, 1^. To see this let us consider a charged particle in the presence of a strong uniform 
magnetic field. The Lagrangian is given by 



L = — X + A-X-V{Xi,X2) (1) 

where V{xi,X2) is an abritrary confining potential. For convenience we shall consider Aq to 
be a harmonic oscillator potential of the form 

V = \ujt (2) 

The energy eigestates of this system lie on Landau levels and at the limit m — > (or 
equivalently strong 5) the system is confined to the LLL. In this case the kinetic energy 
term in (|l]) is negligible. We further use the radial gauge Ai = ^BeijXj. The canonical 
momentum of Xi is then pi = Bx2 and result 

[xi, X2] = [xi, 1^] = ^ (3) 



In the absence of a confining potential all the eigenstates at each Landau level are degener- 
ate. The confining potential lifts the degeneracy. In the lowest Landau level the Hamiltonian 
is 

H=^ia^a + ^) (4) 



where a = + 1x2). The single particle eigenstates are 

n) = ^ at"|0) (5) 
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The corresponding normalized wavefunctions are 



where z = \f^{xi + 1x2). These can be thought of as the coherent state representation of 
(I)- 

Consider now the case of fermions. Since each state can be occupied by at most one 
fermion, the presence of the confining potential selects a unique ground state which is the 
minimum angular momentum state. This is the v = 1 ground state wavefunction 



^!i{xi,...,xn) = n(^i -%)e-2 2...il->l (7) 

i<j 



This corresponds to an incompressible circular droplet configuration of uniform density p 

B_ 

2lT 



This incompressibility is crucial in explaining the experimentally observed gap for u = 1. 



The existence of a similar gap is less obvious in the case of noninteger filling fractions, 
where the LLL is only partially filled. In this case, it is believed that the repulsive Coulomb 
interactions among electrons are important in generating strong correlations which eventually 
produce a new ground state with a gap. Laughlin proposed that such a state, in the case 
V = 2f^y is described by a wavefunction 

^2p+i{x,, ...,xn) = f[{z^ - z,)'^+'e-iS"i l^^l' (8) 

i<j 

Using the connection to the one-component, two-dimensional plasma, Laughlin showed that 
this corresponds to an incompressible droplet of density p = 2n{2p+i) ■ Although not exact 
eigenfunctions, (^) are quite close to the true solutions, at least numerically. They vanish 
quite rapidly if any two particles approach each other, and this helps minimize the expec- 
tation value of the Coulomb energy. Their success lies very much on their short distance 
behavior. 

Our aim is to develop an appropriate coherent state representation for the states of 
the corresponding matrix model and compare them directly to the Laughlin wavefunctions. 
Before we explain this in more detail, we shall give a brief review of the matrix Chern-Simons 
model proposed by Polychronakos. 



3. Chern-Simons matrix model 

The action describing the matrix Chern-Simons model is given by 

S = J dt^Tr{e,h{Xa + i[Ao, Xa])X, + 2Mo - uX'J + - A^^) (9) 

where X^, a = 1, 2 are N x N matrices and is a complex iV- vector that transforms in the 
fundamental of the gauge group U{N), 

Xa^UXaW , ^-^[/^ (10) 

The Aq equation of motion implies the constraint 

G = -iB[Xi, X2] + ^^"^ -Be = (11) 
The trace of this equation gives 

= NB9 (12) 

It is interesting to note that in the absence of \l/'s, the parameter 6 has to be zero. In this case 
the action (^) is that of a one- dimensional hermitian matrix model in a harmonic potential, 
which is known to be equivalent to N one-dimensional fermions in a confining potential 0. 
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Upon quantization the matrix elements of Xa and the components of \1/ become operators, 
obeying the following commutation relations 

1 

{X2)ki\ = -^Sii 6jk (13) 

The Hamiltonian is 

^ = ^(^ + E4a-^) (14) 

where A = ^J^{Xi + iX2). The system contains A^(A^ + 1) oscillators coupled by the 
constraint (p!l|). As explained in 0], upon quantization, the operator G becomes the generator 
of unitary rotations of both Xa and \E'. The trace part (^) demands that NB6 being the 
number operator for \I^'s is quantized to an integer. The traceless part of the constraint 
demands the physical states to be singlets of SU{N). 

Since the Ajj transform in the adjoint and the \I^J transform in the fundamental repre- 
sentation of SU{N), a purely group theoretical argument implies that a physical state being 
a singlet has to contain Nn \l/^'s, where n is an integer. This leads to the quantization of 
B6 = k. We shall see later that in identifying this matrix model with a fermionic system, k 
has to be an even integer. 

Explicit expressions for the states were written down in |^ . (Essentially equivalent results 
were also obtained by Polychronakos [^].) The ground state being an SU{N) singlet with 
the lowest number of A^'s is of the form 

= [e^-^-vl/J^(vl/tAt),,...(vI/tAt^-i),^]>) (15) 
where |0) is annihilated by A's and ^^'s, while the excited states can be written as 

N-l 

|v^,.,) = n (TrAty'[e^-^-vl>t^(v[/tAt),^...(vI/t^tA^-i),^]>) (16) 

i=l 



4. Coherent state representations 

We now present a coherent state representation for the matrix states ( p!5| , |T6|). In doing 
this one has to make a choice of coordinates, which eventually parametrize the phase space 
of the underlying one-dimensional system. 

From the matrix model point of view, a natural choice would be to diagonalize A and 
interpret its eigenvalues as phase space coordinates of particles. We refer to this as the 
^-representation. This is essentially a generalized complex random matrix model. 
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Another choice would be to diagonahze Xi and interpret its eigenvalues as one- dimensional 
coordinates. This is the X-representation. In this, the elements of X2 are canonically conju- 
gate to Xi. Once the X-representation of the matrix states ( pTS] , |I6[) is derived, we can then 
use the usual coherent state representation to express the wavefunctions in terms of phase 
space coordinates. 

Clearly in the two cases. A- and X-representation, the notion of phase space coordinates 
is different; as a result we derive different expressions for the corresponding wavefunctions. 

Below we present in detail the two representations and compare the results to the Laughlin 
wavefunctions. 

A — representation 

We define the coherent state | Z, 0) such that 

Amn\Z,(f)) = Zmn\Z,(f)) 

= 0„ \Z,(f)) (17) 

where Z is a complex N x N matrix and is a complex vector. Let us consider the matrix 
ground state of the form 

|2p+ 1) = [e^-'-vI/t^(vI/tAt),,...(vl/tAt^-i),^]2p|o) (18) 

and reexpress its scalar product in terms of the coherent state wavefunctions using the 
completeness relation 

N 

(2p+ l|2p+ 1) = J \{dZi^ dZ*^ dcPi #;(2p+ l|Z,0)(Z,0|2p+ 1) (19) 

Before we talk about the general case, it is interesting to demonstrate how the coherent state 
representation works for p = 0. The scalar product ([T9|) can be written as, 

. N 

(1|1) = / l[dZ,jdZ*je-^''^^^ / d<Pid4)ie-~^'*' (20) 



Since the (p integration is trivial, eq. (poD defines essentially a complex random matrix model 
with probability distribution e"^'''^^^. A detailed analysis of this model has been given in 



Z is a complex matrix that can be diagonalized as, (nondiagonalizable matrices form a 
set of measure zero) 

Z = XEX-^ (21) 

where E is diagonal with En = Zi. Integration over the nondiagonal part of Z in (pO]) gives 
the following result |l^ 

(1|1) = j l[dzdz^e-^''^"l[\z,-z,\' (22) 

i=l i<j 



up to normalization factors. One can immediatelly recognize the integrand of (P^ as 
the probabihty distribution corresponding to the u = 1 Laughhn wavefunction. In A- 
representation we have identified the phase space coordinates of the fermions with the eigen- 
values of the matrix Z in (ITT 



We would like now to extend this approach in the case p 7^ in the presence of the extra 
\l/ degrees of freedom. This can be viewed as a generalized complex random matrix model. 
(An attempt for a random matrix formulation of the Laughlin theory of the fractional QHE 
has been proposed by Callaway As shown in [10| any complex regular N x N matrix 



X can be expressed in one and only one way as 

X = UYV (23) 

where U is a unitary matrix, y is a triangular matrix with Yij = for i > j and Yu = 1 
and is a diagonal matrix with real positive diagonal elements. Using this particular 
parametrization we find 

< 2p+l\Z,^>= [6^-^-0,,(Z0),,(Z^0),,...(Z^-V)^.]'^e-^^^'^e-i^^ 

= [e''---'''^i,iUYEY'^U~^(P)i,{UYE^Y'^U~'(P)^,...iUYE^-'Y-^U~^^)i^]^^Pe~^^^^^ 

= [det([/y)]2p Y[{zi - Zj^P n ^■''e~2^'^'^e-2'^'^ (24) 

i<j i=l 

where ^ = {UY)-^(j). 

Further using (^TD and (|^) we have 



TrZ^Z = Tr{EHEH-^) (25) 

where H = Y^Y and det H = 1 (det Y = 1). 

In performing the integration over 0's it is convenient to change variables from (j) to ^. 
Since U is unitary and detF = 1 

N _ _ N _ 

n d^i d<Pi e-^'^ = n dii dil e-«^« (26) 
1=1 1=1 



Further following one can show that 



dZijdZ*j = Yl dzidzi H I -^^ ~ -^i I ^ 11 

i i<j i^j 

= dfi{z,H)ll\zi-z,\^ (27) 



upto normalization factors. 
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Putting everything together we find 



r _ N _ 

{2p +l\2p+l)=J df^iz, H) n - ;^,f ^+^e-^(^™-^)de,cie, (28) 

Integration over ^ and H would produce a quantity that depends only on Zi, Zi. It is 
clear from the above expression that in the absence of the ^ integration, the integration 
over the nondiagonal elements of Z would produce a probability distribution similar to the 
corresponding Laughlin one. Since however ^'s couple to H and H in turn couples to Zi and 
Zi, this integration will not necessarily produce a probability distribution which agrees with 
the Laughlin one. The integration over and if is a lot more involved now and it is hard 
to extract a closed expression for arbitrary A^. However, as we shall see one can find general 
features of the probability distribution, which do not agree with what one would expect from 
the Laughlin distribution. 

We first perform the ^-integration by introducing a source term as follows 

_ N 

d^^dl,X{{Ui)^^e-^''^ 

i=l 



|J,J=0 



^99 



Eq. (|28D can now be written as 

i<j i=l ^-^i ^"^i " i^j 



(2p+l|2p+l) = / X{d^^d-z^X{\z,-z,\'^^'X{{-^^f' j nrfi^.,e-^(™-^)+'^"^-^^ 



J,J=0 

(30) 



We are now left with the H integration 

/ = / n rfi^i,e-^(^^^^"')e^"^"'^ (31) 

To do this we follow an iterative procedure as in [0. At each step of the iteration we 
integrate over the variables of the last row and column of H and thus decrease by one the 
size of the matrix. In the absence of the source terms, the structure of the reduced matrix 
remains the same and this produces a simple recursion formula, which eventually leads to 
(^). This is not quite the case when there are source terms, as a result there is no simple 
recursion formula. 

The iteration procedure is defined as follows: Let H' , E', etc., be the relevant matrices of 
order n and H, E, be those obtained from H', E' by removing the last row and last column. 
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Greek (Latin) indices run from 1 to n (n — 1). Let A'^,^ be the cofactor of H'^^ in H' and Ajj, 
the cofactor of Hij in H. Let gi = H[^. Because det H' = det H = 1, the following relations 



are true 10 



A' — w'-i 



A^ = i7;„A,,-^^:^,A 



* „ A. Ik 
ij 



Lk 



Ag = A,,Aik - AikAij 



Let 



Using (|3^) we find 



= Tr{E'H'E'H -') - J'R -'j' = <l>° - J'H "V 

-{g\H-'jJH-'\g) + g*{H-^J)iJN + JN{JH-%gi 



Integration over (7's produces the following result 



l<ij^j<N~l 



det(Xy - JiJj) 



where 



X,, = {H-\E - zm)H{E - zm)H~^),^ - Xjj = - J J, 
The expression ( |35D can be further simplified using the following relations 



Further 



det(X° - J J) = eTrl°g(^°--^'^) = gTr[logXO+log{l-(XO)-ijJ)] 

det(X°) ei°^^"(^°)-^^ =E\z^- zn\' E MX^^^Jn 

i<N 



l<Lm<N-l 



where we have used that 



detX° = det{H~\E - zn)H{E - zn)H-^) = [] 

i<N 



Zi - Zn\ 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 
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Substituting in 



we find 

P-\ZN? 



n 



where 



i<i,j<N-l 



n 



e 



l<if^j<N-l 



M) 



JnJn 

e i-A^ 



E 



J, 



2;^ - zn) 



J. 



Zn) 



(40) 



(41) 



Already this resuh can be used to evaluate (|28|) for the simple case of a 2 x 2 matrix model. 
Although this is a rather trivial case, it is worth presenting it, since it highlights properties 
of the probability distribution which are not in agreement with the Laughlin one. 

To simplify the calculation let us further choose p = 1. For the N = 2 case after the first 
iteration we find from (|4^) 

12 U,|2 

J2J2 



\zi-Z2\^ l-M 



(42) 



where M = 1 — j^^^^^.The functional diferentiation with respect to J's, eq. (|29D , gives 



' d d 
dJidJi 

-- 2 



d d \ 
dJ2 dJ2 ) 
d d 



I J,J=0 



3 Ji Ji 



\Zl - Z2\ 



Using 



1) and 



\Z\ - Z2\ 

we find 



dJidJiJ 
24 



4 + 



\Zl - Z2\ 



1 - 
+ 



I21-22P 

48 



3 lJi,Ji=0 



1^1 ~ 



(43) 



(3|3) 



dzidzidz2dz2e 



-|-lP-k2| = 



Fl - Z2\ 
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Fl - Z2\ 



\Zl - Z2\ 



(44) 



The first term corresponds to the probability distribution for the u = 1/3 ground state 
Laughlin wavefunction. There are extra terms though, which are dominant at short distances 
as zi ^ Z2. In this simple case we find that the distribution emerging from the matrix model 
has a long distance behavior similar to the corresponding Laughlin one, but its short distance 
behavior is quite different. We shall now argue that this behavior prevails for any A^. 



It is clear from (^0]) that the first step of the iteration produces an expression which is 
not quite similar to the original one if J's 7^ 0. As a result the integration over H'^^ at 
each subsequent step of the iteration becomes quite involved. Although it is very hard to 
derive an exact expression for / as a function of Zj's, it is quite straightforward to explore 
its dependence on \zfyf_i — zi^\. This is sufficient for example to get information about the 
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short and long distance behavior of the probabihty distribution as \zn_i — Z]\[\ << 1 and 
\zn-i — Z]\f \ >> 1 respectively. 

The dependence of I on \zn-i — zn\ comes first from the overall factor 1/ Yli^N ~ ^nI"^ 
and second from the factor M, see eq. (ffOl), (pj). Expanding M we find 



Jn-iJn-1 . Jn-I Jj Ji Tj-l Jn-1 



|2;Ar_i - 2;Ar|^ (ZAr_l - ^A?) [Zj - Zn) {Zi - Zn) {Zn-1 - Zn) 

+ E , , (45) 

It is clear now that in order to explicitly demonstrate the I^at — -2iv_i| dependence we need 
to evaluate the functional derivatives with respect to and Jn-i, as in (^0|) . We find 

TT f d 9 \' , e-Sj--P 5 B* 



dJidJiJ ' nj<Ar ki — ^^A^P Zat-i — Zat Zn-1 — Zn 

C D D* 

+ 7 + 



\ZN-1—Zn\^ {ZN-l—ZNy {zn-I — ZnY 

I ( E E \ F , 

\Zn.I- ZnY \zn~i- Zn Zn-I-Zn) \ZN-l-ZNr 

where A, B, C, D, E, F are in general functions of z's containing factors of the form 
{zi — Zn), {zi — zn-i), {zi — Zj) and their complex conjugates, where i < N — 2. Substituting 
this expression in (jS^) we find a probability distribution which is dominated by a term 
proportional to \zn — ZN-if when \zn — zn-i\ » 1, similar to the z/ = 1/3 Laughlin 
distribution. However when \zn — zn-i\ << 1 the matrix distribution is dominated by a 
term proportonal to \zn — Zn-i\'^- The long and short distance behavior for any is the 
same as the one found in the simple N = 2 case, eq.(^). 

It is straightforward now to see that similar results can be derived for any p 7^ 0. In 
particular the long distance behavior of the probability distribution is that of the v = 
Laughhn distribution while the short distance behavior is that of a. u = 1 Laughlin 
distribution. 

Further the probability distribution cannot be factorized as \E'*\E' where \E' is the corre- 
sponding many-particle wavefunction, which is both antisymmetric and holomorphic in z's. 
This indicates that the identification of the eigenvalues of the matrix Z with the actual 
holomorphic coordinates of fermions may not be appropriate. 

X — representation 

Going to an X-representation first we derive one-dimensional fermionic wavefunctions, 
where the coordinates have been identified with the eigenvalues of the matrix Xi, then 
transform the wavefunctions to the coherent state representation. 
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We define the state \X, 0) such that 
Xi|X,0) = X\X,(f: 



IX 



Using (|47| ) we find 



(47) 



(48) 



where 



(O|X,0) 
A, 



2^ BdX,/ 



(49) 



Since (|^) is completely antisymmetric in the i^-indices, the differential operator P^O" 
duces a nonzero contribution only if it acts on the ground state wavefunction (0|X, 0). We 
then find 



{2p+l\X, 



N-li\ -\2p -Tt^BX 



(50) 



X being a hermitian matrix, it can be diagonalized by a unitary transformation 

X = UxU~^ Xij = Xi6ij (51) 

Using this in (|50D we find 



(2p+l|X, 



N-1tt-1i\ ]2p TrifiX 



e ^"2^^ e 2^ 



\2p 



N 



(52) 



i<j 



i=l 



Using (|5^ ) we can express the scalar product of the \2p + 1) state as 
{2p+l\2p+l) = j[dXi^]d^id(t)i{2p+l\X,(t)){X,(t)\2p+l) 



[dXij] W{xi - Xj) 



*P-TtBX^ 



n2p 0</, 



(53) 



Since U is a unitary matrix the integration over 0's completely decouples (unlike the case in 
eq.(^). Further § 

[dXij] = dxi \{{xi - Xjf[dU] (54) 
where [dU] is the Haar measure. Integration over the nondiagonal elements of X gives 



{2p + l\2p + 1) ~ / dxi Wix, - a;,)^P+2g-B^^x: 



(55) 



i<3 

12 



This is the probabihty distribution for the one-dimensional Calogero ground state wavefunc- 
tion Yli^j{xi — Xj)'^P~^^e~^^i^^^^ . This is not surprising; it was aheady indicated in [Q, |12 
that the Chern-Simons matrix model is equivalent to the Calogero model. 

Let us now use a coherent state representation for the Calogero wavefunction (xi, XAr|\E') 
Yli^j{xi — Xj)'^P^^e~^^i^^^'^. The coherent state representation of any wavefunction (x|\E') 
can be written as 

= J dx{z\x){x\'^) (56) 
where z\z) = z\z) and z = J ^{x + iy)^ [x, y] = Using 



z\x) = e-f-'e^^"^e-^"/2e"l^l'/2 (57) 



we find 



= / dxe'^^"""^^' fix) (58) 



where /(x) is given by (x|\E') = f{x)e ^^^1"^, In evaluating (^) we expand f{x) around 
z/^/2B. 



fc=0 



e 2 2^ 



-e-'-^[e^hf(x)l^,^^ (59) 



The coherent state representation of the many-body Calogero wavefunction is thus 



Zn 



For p = (0) gives the v = 1 Laughlin state. For p 7^ the coherent representation of 
the Calogero state has a long distance behavior similar to the v = Laughlin state, but 
a different short distance behavior. This connection between the one-dimensional Calogero 
wavefunction and the Laughlin states has already been noted in [13|. 



5. Summary and discussion 

In an attempt to clarify the exact correspondence between the Chern-Simons matrix 
model introduced by Polychronakos and the fractional QHE at filling fraction u = 1/m, 
as described by Laughlin wavef unctions, we have derived the matrix model wavefunctions 
using a coherent state representation. We have presented two different coherent state repre- 
sentations, each one implementing a different choice for the phase space coordinates of the 
underlying one- dimensional fermionic system. 
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In the ^-representation, the eigenvalues of the matrix A are identified with the phase 
space coordinates z of the fermions, while in the X-representation, the eigenvalues of the 
matrix Xi are identified with the one-dimensional coordinates x of the fermions. 

Both choices give identical results when p = 0, or equivalently 6' = in (P). The 
corresponding wavefunction is identical to the u = 1 Laughlin wavefunction. 

For p the two representations give different results. Although the explicit expressions 
for the probability distributions are different, they share some common features. They both 
have the same long and short distance behavior. Comparing them to the corresponding 
u = Laughlin distributions, we find that it is only the long distance behavior which is 
in agreement. The short distance behavior does not agree with the Laughlin one. 

A noticeable difference between the two representations is that the X-representation leads 
to a holomorphic wavefunction with antisymmetric properties, while this is not possible for 
the ^-representation. 

As we mentioned earlier there is an ambiguity in introducing electron coordinates in 
the matrix model. This has to do with the choice of coherent state representation. In 
this paper we analyzed two particular coherent state representations, which seem to be 
the natural choices from the matrix model point of view. In both cases we find that the 
emerging wavefunctions do not quite agree with the Laughlin one. Although this by itself 
does not prove that the original matrix model is not equivalent to the Laughlin theory for 
the V = 1/m fractional QHE, it makes the precise correspondence between the two models 
less transparent. One can argue that there may exist another coherent state representation, 
corresponding to another choice for the electron coordinates, such that there is agreement 
with the Laughlin wavefunction. However, for the matrix model to be truly useful in the 
context of QHE, this new coordinate choice should be easily identifiable. 
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